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Abstract: Integrating seminal ideas of London, Feynman, Uhlenbeck, Bloch, 

Bardeen, and other illustrious antecessors, this paper continues the development of 
^ an ab initio theory of the A transition in liquid *He. The theory is based upon vari- 

ational determination of a correlated density matrix suitable for description of both 
I normal and superfluid phases, within an approach that extends to finite tempera- 

^ ' tures the very successful correlated wave-functions theory of the ground state and 

CJ ■ elementary excitations at zero temperature. We present the results of a full opti- 

mization of a correlated trial form for the density matrix that includes the effects 
both of temperature-dependent dynamical correlations and of statistical correlations 
corresponding to thermal phonon/roton and quasiparticle/hole excitations - all at 
the level of two-point descriptors. The optimization process involves constrained 
functional minimization of the associated free energy through solution of a set of 
Euler-Lagrange equations, consisting of a generalized paired-phonon equation for the 
structure function, an analogous equation for the Fourier transform of the statistical 
exchange function, and a Feynman equation for the dispersion law of the collective 
excitations. Violation of particle-hole exchange symmetry emerges as an impor- 
tant aspect of the transition, along with broken gauge symmetry. In conjunction 
with a semi-phenomenological study in which renormalized masses are introduced 
for quasiparticle/hole and collective excitations, the results suggest that a quantita- 
tive description of the A transition and associated thermodynamic quantities can be 
achieved once the trial density matrix is modified - notably through the addition of 
three-point descriptors - to include backfiow effects and allow for ab initio treatment 
of important variations in effective masses. 
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1. INTRODUCTION 



The remarkable coherent collective property of superfluidity, observable in liquid 
'^He below the famous A transition, was discovered by Kapitza [1] and Allen and Mis- 
ener [2] in 1938. Almost immediately, Fritz London [3,4] proposed that the transition 
to the supcrfluid phase (Helium II) reflects a process that is fundamentally analogous 
to Bose-Einstein condensation of an ideal Bose gas: 

I recently realized that ... some support could be given to the idea that 
the peculiar phase transition ("\-point"), that liquid helium undergoes at 
2.19 K, very probably has to be regarded as the condensation phenomenon 
of the Bose-Einstein statistics, distorted, of course, by the presence of molec- 
ular forces and by the fact that it manifests itself in the liquid and not in 
the gaseous state. 

Based on this idea, Tisza [5,6] introduced a successful two-fluid model of the macro- 
scopic behavior of Helium II; a similar model was developed independently by Landau 
[7] on different grounds. 

Here we shall report continued progress toward a quantitative microscopic un- 
derstanding of the A transition in liquid ^He in terms of Bose-Einstein condensation 
of strongly interacting bosons. Carried out within the semi-analytic framework of 
correlated density-matrix theory, this effort aims at a concrete atomistic realization 
of London's hypothesis as further shaped by Feynman [8] . In developing the theory, 
we have encountered a number of signposts erected by great physicists, suggesting 
that the path being followed is the right one. As a preface to the presentation of new 
results, it will be instructive to re-examine some of these seminal insights. 

Working in conflguration space, the density matrix of a system of N identical 
bosons can be written in the general form 

W(R, R') = /-^$(R)Q(R, R')$(R') , (1) 

where the coherence factor $(R') becomes the ground-state wave function in the 
limit of zero temperature, (5(R, R') is an incoherence factor that is non-separable in 
the conflgurations R = (ri, . . . , r^v) and R' = (r^, . . . , r^) and describes the effects 
of real excitations, and is a normalization constant. It is convenient to split the 
incoherence factor into two parts, (5(R, R') = Qcon(R-,^')Qqp(Ri^')- The factor 
Qcoii embodies the effects of real, thermal collective excitations (phonons and rotons), 
while Qqp accounts for the effects of real quasiparticle excitations (or exchange corre- 
lations). In the trial density matrix currently employed, the factors $, Qcoib and Qqp 
are built from two-point functions. In particular, $(R) and QcoiKR) R-') ^ire taken 
essentially as Jastrow products of appropriate dynamical and collective-thermal cor- 
relation functions [9-11] exp['u(|ri — rj|)/2] and exp7(|ri — rj|), respectively, while 
Qqp(R, R') is expressed in terms of a permanent of appropriate statistical correla- 
tion functions [12] Fccdri — r^|). Ultimately, the functions u, 7, and Fee are to be 
determined by functional minimization of the free energy. A description of backflow 
effects requires the inclusion of three-point functions in the coherence factor $(R) 
and the incoherence factor (5(R, R'). 

Decomposition of the factors $(R), $(R'), (5qp(R, R')? and (5coii(R, R') into 
products of two-point functions (and, as needed, three- and higher-point functions) is 
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motivated by the quantitative successes of the variational wave-function theory of the 
ground state of a strongly interacting Bose system [13]. In this theory, the ground- 
state wave function is written as a Jastrow product of pair correlation functions 
supplemented by a product of triplet correlation functions, and so on (as needed) to 
higher-body products in the Feenberg representation of the exact Bose ground state 
[11]. 

The treatment of (5qp(R, R') in terms of a permanent of two-point functions 
(the statistical correlation functions Fgcdri — ^'jD) exploits and extends the ideas of 
Feynman [8] and the formalism of Ziff, Uhlenbeck, and Kac [12]. For free bosons, 
which is the case studied by Ziff et al., Tcc{r) has a Gaussian form. Feynman used 
this form but inserted an effective mass in place of the bare mass. 

To quote Feynman's abstract. 

It is shown from hist principles that, in spite of the large interatomic forces, 
liquid "^ffe should exhibit a transition analogous to the transition in an ideal 
Bose-Einstein gas. The exact partition function is written as an integral over 
trajectories, using the space-time approach to quantum mechanics. It is next 
argued that the motion of one atom through the others is not opposed by 
a potential barrier because the others may move out of the way. This just 
increases the effective inertia of the moving atom. This permits a transition, 
but of the third order. It is possible that a more complete analysis would 
show the transition implied by the simpliGed partition function is actually 
like the experimental one. 

The strategy now being pursued within correlated density-matrix theory seeks to 
complete Feynman's program using a description that is more flexible yet is rooted in 
the same physical picture in which exchange effects are of paramount importance. As 
the system goes from the normal to the superfluid state, exchange loops (or "rings" 
or "cycles") of macroscopic length come to dominate the behavior of the system; this 
feature is revealed in the development of off-diagonal long-range order of the partial 
distribution function (or correlation function) Gcc{f) that characterizes exchange 
processes. The path-integral Monte Carlo approach of Ceperley and Pollock [14,15] 
is also based on Feynman's ideas on liquid "^He, played out within Feynman's space- 
time approach to quantum mechanics. Quite naturally, a similar picture emerges in 
this work - involving a kind of percolation transition in which exchange paths can 
grow to macroscopic size. 

The entropy corresponding to the chosen trial density matrix is determined us- 
ing a replica ansatz (borrowed from the statistical physics of disordered systems), 
together with the so-called diagonal approximation, which implements a separability 
hypothesis analogous to that applied by Campbell and Feenberg [9] in the paired- 
phonon theory of the ground state. The work of Ziff, Uhlenbeck, and Kac [12] 
provides the pivotal clue for constructing a trial density matrix that is applicable in 
both normal and condensed phases and across the phase transition. 

The presence of a permanent in the trial density matrix introduces complica- 
tions in the evaluation of spatial distribution functions and structure functions and 
of the internal-energy and entropy components of the free energy. However, a pow- 
erful formalism for handling exchange correlations comes to the rescue: the technical 
problems created by the permanent can be solved by appealing to (and extending) 
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the methods of Fermi hypcrncttcd-chain theory invented by Fantoni and Rosati [16] 
and Krotscheck and Ristig [17]. Fermi-hypernetted chain theory teaches us how to 
evaluate correlated expressions involving determinants of two-point statistical func- 
tions, and this knowledge is readily adapted to Bose statistics. Complications also 
arise in the evaluation of the entropy using the replica ansatz, and these can be 
largely overcome with the aid of methods developed by Hiroike [18] for classical fluid 
mixtures. 

The entropy, calculated in diagonal approximation, decomposes into independent 
contributions from phonon-roton excitations and from two kinds of quasiparticle ex- 
citations. The internal energy also decomposes into contributions from these three 
excitation branches, plus a correlation-energy term that becomes the ground-state en- 
ergy at zero temperature. Within the correlated density-matrix theory, as explicated 
using the hypernetted-chain formalism, it is natural to interpret the two particle- 
type excitations as quasi/io/es and qvLa,siparticles. The quasihole and quasiparticle 
dispersion relations are identical in the normal phase, where particle-hole symmetry 
is maintained, but they differ markedly in the condensed phase, where the exchange 
symmetry is broken. The A transition therefore involves two order parameters: 

(i) A condensation strength parameter Bcc which vanishes in the normal phase and 
measures the breaking of gauge symmetry in the condensed phase (also occurring 
in the noninteracting Bose gas) 

(ii) An exchange strength M which also vanishes in the normal phase and measures 
the violation of exchange symmetry below the A point (occurring only if the 
interactions are strong enough) 

On physical grounds, these order parameters must be coupled, the simplest reasonable 
connection being a straight proportionality M = aBcc, where the coupling parameter 
a is temperature dependent [19]. 

The interpretation of the two quasiparticle branches as quasihole and quasiparti- 
cle modes gains support from an idea voiced by Felix Bloch to the effect that - much 
as in a Fermi system - there might exist quasiparticle and quasihole excitations in a 
system of identical bosons. Bloch's idea (unpublished, but aired at a Stanford semi- 
nar) resonates with Pines' phenomenological arguments [20] that strongly interacting 
Bose and Fermi liquids have much more in common than is normally supposed, a 
claim well documented by experimental findings on the properties of their elementary 
excitations. 

The nature of the proposed exchange-symmetry-breaking phenomenon is illu- 
minated by the existence of a striking parallel with the behavior of a diamagnetic 
material in a magnetic field [19]. A particle-hole exchange field may be constructed 
as the sum of differences between the hole energies and the corresponding free single- 
particle energies, weighted by normalized differences between the occupation numbers 
of the two types of quasiparticle excitations. Pursuing the electromagnetic analogy, 
a macroscopic polarization is induced by this field, which is negatively proportional 
to the field with a diamagnetic-susceptibility coefficient that depends on the occu- 
pations of the two quasiparticle branches. This polarization effect acts to screen 
the exchange field in the condensed phase. Indeed, in analogy with the Meissner 
effect, the condensed phase behaves like an ideal diamagnetic as T approaches zero, 
completely expelling the exchange field from the system. The screening leads to a 
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drastic reduction in the particle-hole exchange contribution to the internal energy 
of the system. This exchange property was considered in detail in Ref. [19], where 
the exchange-symmetry breaking was first proposed as an essential aspect of the A 
transition in liquid ^He. 

The diamagnetic analogy evokes a conviction expressed by John Bardeen [21] in 
recounting the development of superconductivity theory: 

It seems to me that most of those who thought long and hard about super- 
conductivity prior to the discovery of the Meissner effect in 1 933 never got 
over an inner feehng that the reaUy fundamental property of a supercon- 
ductor is infinite conductivity or persistent currents, and this colored the 
way they thought about the subject in future years. While an adequate 
theory must explain both aspects, the diamagnetic approach has been the 
most fruitful in indicating the nature of the superconducting state. 

This deep insight sanctions an analogous description of the superfiuid state of liquid 
^He. 

In advance of a complete functional optimization of the trial density matrix, a 
model calculation [19] has provided strong support for the new microscopic approach 
and especially for the importance of exchange-symmetry breaking in the A transition. 
A standard parameterized form was assumed for the dynamical correlation function 
(i.e. for the pseudopotential u{r)), while - as in Feynman's treatment - the statistical 
correlation function TcdT) was determined through an effective-mass parameteriza- 
tion of the quasihole spectrum. The collective (or phonon-roton) incoherence factor 
QcoiKR-j which is not expected to exert a strong influence on the transition, was 
set equal to unity. A simple and plausible model was introduced for the coupling 
parameter a(T) connecting the order parameter M for exchange-symmetry violation 
to the order parameter Bcc for gauge-symmetry breaking. This model requires, as 
empirical input, the ratio of the speeds of second and first sound in the condensed 
phase. 

Within the semi- microscopic/semi-empirical model so delineated, exchange- sym- 
metry breaking was found to reduce the condensation temperature of liquid '^He from 
the value 3.2 K corresponding to Bose-Einstein condensation of a noninteracting Bose 
gas to a value 2.2 K very close to the experimentally measured A transition temper- 
ature. This prediction is based entirely on the extrapolated behavior of the hole 
chemical potential in the normal phase; it is consistent with the behavior of the 
model in the condensed phase but does not depend on the empirical input for the 
sound-speed ratio. The prediction for the specific heat shows an apparently diver- 
gent behavior at 2.2 K and qualitatively - though not quantitatively - reproduces 
the characteristic A shape of the experimental curve. The latter result contrasts 
with that of Feynman's treatment, where an approximate and incomplete analysis 
of the geometry of exchange processes led to a continuous specific heat, as in the 
noninteracting Bose gas. 

In this paper we continue to build upon the correlated density-matrix theory 
of strongly correlated Bose systems, essentially in the form designed and developed 
in Refs. [22-27]. The trial density matrix for liquid '^He and the consequent struc- 
ture of the energy and entropy (the latter evaluated in diagonal approximation) are 
made more explicit in Sec. 2. The Euler-Lagrange equations for the optimal density 
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matrix without backflow, expressed in convenient form (as presented in Sec. 3), are 
solved numerically. The results of the full optimization are examined in detail in 
Sec. 4. (Additional information is furnished by Ref. [28].) Remaining obstacles to a 
full microscopic understanding of the A transition are discussed in the final section. 
Detailed formulas of the theory are relegated to an appendix. 



2. CORRELATED TRIAL DENSITY MATRIX AND FREE ENERGY 



The theory begins with suitable choices for the coherence and incoherence factors 
appearing in Eq. (1): 

(i) A wave function of Jastrow form 



^R) = u E 



$(R)=exp<^-> zx(|r,-r,|) V , (2) 



(ii) An incoherence factor Q = QcoiiQqp that takes account of real collective excita- 
tions (phonons, rotons) through 



^ r 1 1 1 

Qeoii(R,R') = exp J] I 7 (|r. - r^ l) - -7 (|r. ' - H^i ' ^'jD \ ' 



(3) 



and incorporates the effects of quasiparticle excitations via 

gqp(R,R') = ^^/^e^/-Perm[r,,(|r,^r;.|)+S,,z] . (4) 

The contour in Eq. (4) encircles the origin in the counterclockwise sense, and 
the quantity A is an integral involving Tec and the dd and dc components of 
the radial distribution function [19]. The condensation strength Bcc serves as 
the primary order parameter that monitors the breaking of gauge symmetry as 
one enters the Bose-Einstein condensed phase, being restricted by < Bcc < 1- 
In the normal phase where Bcc = 0, the integration in expression (4) may be 
carried out analytically and the incoherence factor (5qp(R, R') reduces simply 
to the permanent of the statistical functions Fccdri — fjl) alone. 
The free-energy functional associated with the trial density matrix specified by 
Eqs. (2)-(4) has been constructed explicitly [19], with the following general results. 
The internal energy of the A'"-boson system at temperature T is the sum of a cor- 
relation energy Ec (which coincides with the expected ground-state energy at zero 
temperature) and separate contributions from collective, quasihole, and quasiparticle 
excitation branches, 

E = Ec + EcoW + Eqh + Eqp . (5) 

The energy addends -Ecoii, -Eqh, -Eqp account, respectively, for the individual renor- 
malized energy contributions ecoii(fc), ^cc{k), and eo{k) = h^k^ /2m from phonons (or 
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rotons), quasiholes, and quasiparticles of all momenta Kk, 

-Ecoii = ^ eco\i{k)n{k) , (6) 

k 

-^qh = ^ XI ^cc{k)ncc{k) , (7) 

k 

Ec,^ = lY.'oik)n,{k). (8) 
k 

In turn, the functions n{k), riccik), and ndk) are the respective thermally aver- 
aged occupation numbers of these three kinds of excitations, at wave vector k. The 
corresponding excitation energies uj{k), Ucc{k), and ujc{k) are defined through the 
Bose distribution function n{k) — [exp (3u!{k) — 1] ^ and through the analogous Bose 
distributions ricdk) and ndk). 

In the diagonal approximation [24,25], the three types of excitations also con- 
tribute additively to the total entropy Sg, 

S, = Si'"''^ + Si^'^ + , (9) 

where 

5(^0") = kBj2{[^ + ^(^)]^^[^ + ^(^)] - ^(^) In n(/c) } , (10) 

k 

— ncc{k) In rice (k)}, (11) 

k 



and Se is obtained from Eq. (11) by replacing ncc{k) with nc{k). 

The energies E'c, ^co\\{k), and ecdk) and the distributions n(/c), ricdk), and 
nc{k) are to be regarded as functionals of the statistical function Tcdr) and the 
pseudopotentials u{r) and 7(r). In the condensed phase, they are also functions of 
the condensation strength Bcc- Explicit expressions for these quantities are provided 
in the Appendix. 

The two excitation branches having Bose distribution functions Ucdk) and ndk) 
and respective energies uJcdk) and ujdk) arc characterized at small wave number k by 
a quadratic dispersion law or /and an energy gap. In the normal phase where Bcc = 
these two branches merge, i.e., Ucdk) — nc{k) and u>cdk) — ^dk) fo^^ ciny /c > 0. 
However, in the condensed phase where Bcc > 0, this degeneracy may be lifted in 
the presence of interactions between the boson constituents. The two branches then 
follow different dispersion laws. In the spirit of Felix Bloch's conjecture as emmci- 
ated in the introduction, we view the cc and c branches as involving quasihole and 
quasiparticle excitations, respectively [19]. According to this Bloch interpretation, 
the symmetry between particle and hole excitations is broken in the condensed phase 
of a strongly interacting Bose system. Further, as sketched in the introduction, the 
exchange-symmetry-breaking phenomenon can be given a formal expression by de- 
veloping an analogy with the behavior of a diamagnetic material in a magnetic field. 
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The counterpart of the Meissner effect (expulsion of the exchange field) has been 
studied numerically in Ref. [19] within a model in which simple parameterized forms 
are assumed for the dynamical correlations and for the hole spectrum that determines 
the statistical correlations. 

To proceed beyond the model of Ref. [19] toward a systematic microscopic treat- 
ment, we now turn to the problem of full functional optimization of the Helmholtz 

free energy. Due to the formal complexity of this problem, the present implementa- 
tion is limited to the choice (2)-(4) for the trial density matrix, implying the omission 
of three- and higher-body correlation functions and hence neglect of what are tradi- 
tionally called "backflow" effects. 



3. FUNCTIONAL OPTIMIZATION OF THE FREE ENERGY 

The Delbriick-Gibbs-Moliere minimum principle for the Helmholtz free energy 
[23] is employed to determine the optimal density matrix (1) within the class of 
trial functions defined by Eqs. (2)-(4). Variation of the functional F = E — TSe 
constructed from the results (5)-(ll) is, however, constrained by a particle-number 
(or particle-conservation) sum rule. This sum rule takes the form [25] 

7V = SccA^o + J]nec(A;) (12) 

k 

with 

No = N + Y.^,,{k) [Sdc{k) + S';^Xk)\ (13) 

k 

in the condensed phase where < B^c < 1; it of course lacks the NgBcc term in 
the normal phase. The functions Sdc{k) and Sic\k) are components of the structure 
function S{k) and the exchange function Scc{k), respectively [19]. Coupling the con- 
dition (12) to the free energy expression with a Lagrange parameter A, the modified 
functional Fx is varied with respect to the input functions u{r), Tcdr), and 7(r) (or 
their Fourier transforms u{k), Tcc{k), and 7 (A;))- The Euler-Lagrange equations 

=0, #^ = 0, #,=0 (14) 



5u{k) ' 5Tcc{k) 5^{k) 

then determine the optimal pseudopotential, quasihole/quasiparticle spectra, and 
phonon/roton spectrum. 

Explicit expressions for the variational derivatives in Eqs. (14) (again, in the 
absence of backfiow) may be found in Ref. [28] . The results can be recast in the form 
of a generalized paired-phonon equation [9,10] 

S{k) + ^eoik)[S{k)-l]=D{k) (15) 

for the radial distribution function or the associated static structure function S{k), 
an analogous equation 

Sccik) + ^eo{k)Scc{k) = D,,{k) (16) 
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for the statistical exchange function Gcc{f) or the corresponding Fourier transform 
Scc{k), and a Feynman equation [29,22] 

(17) 

for the coUective excitations with the dispersion law uj = uj{k) at temperature T = 

{kB0)~^- Equations (15) and (16) involve the generalized structure functions S{k) 

and Scc{k) along with the functions D{k) and Dcc{k). The Appendix and Ref. [28] 
furnish detailed formulas and relations for these quantities. 

To solve the paired- phonon equation (15), one customarily expresses it - equiv- 
alently - as a generalized Bogoliubov equation 

iv\k) = eo{k) [eo{k) + 2vp-h{k) - 2<„ii(^)] , (18) 

where the collective potential is given by fcoii(^) = —2eo{k)n{k){n{k) + l)/S^{k) and 
the particle-hole potential Vp-h{k) is taken from Ref. [28]. 

It is instructive to note that the correlated density-matrix theory may be fruit- 
fully viewed as a pair-density energy-functional approach [30]. The theory provides 
for systematic construction of a universal energy functional E = E[G{r),Gcc{r), Bcc, p] 
which depends on one-body quantities (particle density p and condensation strength 
Bcc) o,nd on two-body densities (the radial distribution function G{r) and the sta- 
tistical exchange function Gcc{r)). Similarly, the internal energy, the entropy, and 
the Helmholtz free energy are all functionals of these one- and two-body quantities. 
Correlated density-matrix theory may be therefore interpreted as a generalization of 
the Kohn-Sham approach [31] that explicitly acknowledges the role played by the 
dynamical and statistical correlations present in the many-body system. 

We further observe that the generalized paired-phonon equation (15) may be 
transformed into a non-linear Schrodinger equation [28,30]. This equation has the 
form [28] 

S^-^A + v{r)+w{r)+Wcoi\{r)+w^p{r)^ = 0, (19) 

where v{r) is the bare ^He-^He interaction and the induced potentials w{r), Wcou{r), 
and Wqp{r) are functionals of G{r) and Gcc{f) and functions of particle density p and 
temperature T. 



2^{k) 



4. NUMERICAL RESULTS AND DISCUSSION 

Evaluation of the various components of the structure functions S{k) and S{k) 

• 

and the statistical functions Scc{k) and Scc{k) is carried out with the hypernetted- 
chain techniques [25,32], neglecting elementary diagrams. The Bogoliubov equation 
(18) is employed to solve the Euler-Lagrange equations efficiently by means of a suit- 
able iteration procedure. In the normal phase, the sum rule (12) serves to determine 
the chemical potential p,cc = —^cciS^) of the quasiholes (and quasiparticles) . The con- 
densation strength characterizing the condensed phase is obtained from relation (12) 
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in conjunction with the assumed optimization condition dFx/dBcc = 0. It must be 
noted, however, that the latter prescription is not entirely adequate; in future work 
it should be replaced by a more appropriate renormalized Hartree equation. The 
HFDHE2 Aziz potential [33] is adopted for the interaction between the '^He atoms. 
The ^He atomic bosons are confined in a cubic box at a uniform density given by the 
experimental saturation density p{T) of the liquid at temperature T, starting with the 

value p = 0.02185 A"' at zero temperature. Reference [28] contains further informa- 
tion on the numerical procedures applied. Optimal versions of the various quantities 
have been calculated in the temperature range T — [0 K, 5.5 K]. Figures 1-7 display 
a selection of the numerical results. 
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Figure 1. The optimal structure function S{k) at various temperatures and at the 
experimental saturation density. 



Figure 1 presents the results for the structure function at several temperatures. 
As expected, the calculated data look very similar to earlier results [23] , demonstrat- 
ing that the exchange correlations (while certainly important for other quantities) 
affect the spatial structure of the ^He liquid only in its details. The variation of 
the low-A; behavior of the structure function reflects the decrease of the velocity of 
first sound with increasing temperature, as has been discussed in detail in Ref. [23] . 

The maximum of S{k) at about /c ~ 2 A shows a monotonic decrease with in- 
creasing temperature. This predicted behavior is in accord with experimental mea- 
surements in the normal phase, but it conflicts with experiment in the condensed 
phase, where the maximum is found to increase with temperature. The generalized 
Feynman relation (17) yields theoretical results for the excitation energies of the col- 
lective phonon/roton branch. The optimal data plotted in Fig. 2 shows the familiar 
dependence of uj{k) on wave number k: the linear dispersion law for phonons at 
low momenta and the typical roton minimum at A: ~ 2 A. However, since backflow 
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contributions have been sacrificed in favor of a closer examination of the effects of 
exchange correlations, the roton minimum is more than twice what is seen experi- 
mentally. Consequently, the number of elementary excitations increases only slowly 
with increasing temperature. To examine more closely the implications of this be- 
havior, let us assume for the moment that Feynman's treatment [34] is applicable 
and that the superfluid density ps of the ^He liquid is approximately described by 
the expression 

Ps = ^-'^Y.^o{k)[l + n{k)]n{k) . (20) 

k 

If the distribution n{k) takes small values, then the excitation energy u{k) is large 
and the sum in Eq. (20) is small; therefore the superfluid density is close to unity. 
Accordingly, high temperatures are needed to make Ps become small. This rough line 
of reasoning already indicates that the current realization of density-matrix theory 
will predict a value for the normal-to-supcrfiuid transition temperature (i.e., the 
temperature T\ of the A transition) that is significantly higher than the experimental 
value. 

This feature of the calculation is clearly evident in the results for the respective 
excitation spectra uJcc{k) and ujc{k) of the quasiholes and the quasiparticles. Figure 3 

shows their dispersion curves for particle density p = 0.02185 A ^ at three differ- 
ent temperatures. The excitation energies are seen to exhibit a minimum around 
the roton minimum for temperatures below T < 3.1 K. In contrast to the case of 
quasiholes, the quasiparticle spectrum u!c{k) is restricted to wave numbers k > ko, 
excitations being forbidden at smaller wave numbers. The threshold ko for the onset 
of quasiparticle excitations decreases slowly with T at very low T, but it decreases 
very rapidly for T > 2.8 K and vanishes at T ~ 3.4 K. At higher temperatures, the 
quasiparticle and quasihole dispersion curves coincide and show a monotonic increase 
with increasing wave number k. 

The low-temperature behavior of the excitation energies can be extracted ana- 
lytically from the Euler-Lagrange equations. One readily finds 

a;cc(A;) = 2[eo(A;)-A] , k<ko, (21) 

uj,{k) = u,,{k) = [eo{k) - X]S-\k) , k>ko, (22) 

where the threshold ko is determined by 2S{ko) = 1. The chemical potential of the 
quasiholes is therefore ficc = ~<-^cc{0) = 2A in the limit T — 0. 

The numerical results for the optimal potential pcc in the temperature range 
< T < 5.5 K are plotted in Fig. 4. As the temperature is lowered from the "large"- 
T side, ficc{T) < increases until it vanishes at a predicted transition temperature 
Tx ~ 3.4 K. In the temperature range T >T\, the liquid is in the normal phase where 
Bcc = 0. For Bcc > 0, there is effectively a potential barrier for quasiholes, formed by 
the sudden decrease and subsequent slow rise of Pcc as the temperature is decreased 
from T < 3 K. As T goes to zero, the hole chemical potential vanishes again and 
we have /Ucc(O) = 2A = 0. However, this low-T behavior of ^cc is an artifact of the 
variational condition dF\/dBcc = that we have imposed on the Bcc parameter. We 
expect that an improved Hartree condition on the condensation strength will produce 
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Figure 2. The optimal phonon/roton excitation energy uj{k) as a function of wave 
number k at various temperatures and at the experimental saturation density. 
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Figure 3. The optimal quasihole energy U3cc{k) and quasiparticle energy U3c{k) as 
functions of wave number k at various temperatures and at the experimental satura- 
tion density. 
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a quasihole chemical potential that increases nionotonically with temperature in the 
condensed phase and vanishes at the A transition point. 

The numerical predictions for the excitation energies uJcdk) and LVc{k) themselves 
are very large, implying that the associated Bose distributions ncc{k) and nc{k) are 
unrealistically small. The particle- number sum rule (12) then leads inevitably to 
a large critical temperature Tx — 3.4 K, even larger than for a Bose gas of non- 
interacting ^He atoms and much too high compared with the experimental A-point, 
Tx — 2.18 K. It would seem that the present ansatz for the density matrix, devoid 
of backflow, is not flexible enough to produce a sufficiently large effective mass of 
the quasiholes in the medium and a concomitant reduction of the excitation energy 
Ucc{k). This perspective is in harmony with Feynman's argument that the motion 
of one helium atom through the others in the liquid increases the effective inertia 
of the moving atom [8]. The parameterized model of Ref. [19] also lends support to 
this assessment: it involves a large quasihole effective mass m* ~ 1.5m and yields a 
theoretical transition temperature Tx — 2.2 K, in quantitative agreement with results 
derived from the present formulation by simulating the effects of backflow (see below). 

Figure 5 plots the calculated optimal data for the internal energy per particle 
E/N as a, function of temperature at the fixed density p = 0.02185 A and at the (T- 
dependent) experimental saturation density. In the temperature range < T < 2 K, 
the theoretical predictions (at saturation density) lie higher than the experimen- 
tal results by about 2.5 K. Numerical calculations within correlated-basis-function 
theory [36,37] as well as variational Monte-Carlo studies [38] at zero temperature 
strongly suggest that the energy excess is due in large part to the omission of triplet 
pseudopotentials u{ri,rj,rk) from ansatz (2). At zero temperature, these terms es- 
sentially account for the energetic effect of backfiow. As expected, the high energy 
per particle obtained in the present treatment is accompanied by a large theoretical 
value for the transition temperature, Tx ^ 3.4 K. 

The Bose-Einstein-condensed phase and the associated violation of gauge sym- 
metry may be characterized by the primary order parameter Bcc, which we have 
termed the condensation strength, or alternatively by the non-zero condensate frac- 
tion no{T). In the temperature interval < T < 3 K, the B^c order parameter 
depends rather weakly on T, decreasing linearly from unity at T = to about 0.8 
at T = 3 K. For slightly higher temperatures, a drastic reduction sets in and Bcc 
vanishes at Tx — 3.4 K. In this narrow region of rapid falloff of Bcc, the behavior of 
the solutions of the Euler-Lagrange equations depends sensitively on the properties 
of the quasihole and quasiparticle excitation branches, which tends to destabilize the 
iteration process we have employed (for details, see Ref. [28]). 

The condensate fraction no{T) has been calculated within the hypernetted-chain 
formalism of Ref. [39]. The numerical results are plotted in Fig. 6, where they are 
compared with the path-integral Monte Carlo results of Ceperley and Pollock [14]. 
At zero temperature, we reproduce the value no(0) ~ 9.8% obtained in an earlier 
application of correlated density-matrix theory and reported in Ref. [40] . This value 
is satisfactorily close to the estimate given by the Green's function Monte Carlo 
approach [41], ?io(0) ~ 9.0 ± 0.3%. With increasing temperature below 2 K, the 
optimal results for no{T) show a gradual linear decrease similar to that suggested 
by the path-integral Monte-Carlo data. However, the expected precipitous decline 
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T[K] 

Figure 4. The optimized quasihole chemical potential //cc versus temperature at 
density p = 0.02185 A (diamonds) and at the experimental saturation density 
(crosses). 
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Figure 5. Optimal theoretical results for the energy per particle as a function of 

temperature T at p = 0.02185 A (diamonds) and at saturation density (crosses). 
The squares mark experimental data taken from Ref. [35] . 
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T[K] 

Figure 6. Optimal theoretical results for the condensate fraction no versus temper- 
ature at density p = 0.02185 A (diamonds) and at experimental saturation density 
(crosses). Shown for comparison are data from path-integral Monte Carlo [14] and 
Green's function Monte Carlo [41] studies (denoted respectively by squares with error 
bars and by the asterisk). 

toward a vanishing condensate fraction is not seen until a rather large temperature is 
reached, T ~ 3.2 K. This behavior, again pointing to a normal-superfluid transition 
temperature that is too high, is consistent with the above findings for the excitation 
spectra, internal energy, and condensation strength. The results of the current treat- 
ment for other observable quantities, including the entropy and velocity of sound, are 
discussed in Ref. [28] . 

To incorporate backflow effects into the present framework, we must generalize 
the formulae (2)-(4) defining the choice of trial density matrix. This may be done by 
supplementing the pseudopotentials u{r) and 7(r) by triplet (or three-point) factors 
and by generalizing the permanent in expression (4) . Alternatively, we may introduce 
higher-order correlations through the device of shadow wave functions [42,43]. Either 
option is both formally and computationally demanding, especially since elementary- 
diagram contributions should be included for consistency. In advance of such a pro- 
gram, it is useful to adopt a more phcnomenological approach: we shall assume that 
"backfiow" has the principal effect of renormalizing the boson mass m appearing in 
the collective and quasiparticle and quasihole terms of the energy (5). Restricting 
ourselves to a treatment of the normal phase, we accordingly modify the functionals 
(6)-(8), derived from first principles, to read 

£^coii = V -^econ{k)n{k) , (23) 

k ""coll 
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Figure 7. The chemical potential jicciT) of the quasiparticles (quasiholes) in the 

normal phase as function of temperature T at density p = 0.02 185 A for different 
values of the effective mass m* 
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(A;) + eo(A;)]ncc(A;) 



(24) 



In general, the effective masses m*^^^ and m*p depend on temperature, density, and 
wave number and differ from each other. 

As anticipated, a revised optimization based on the modified functionals (23) and 
(24) does indeed reveal a strong sensitivity of the theoretical transition temperature 
T\ to variations in the effective mass. Figure 7 collects numerical results for the 
optimal chemical potential of quasiparticlcs/quasiholes in the normal phase at 

density p — 0.02 185 A , for a selection of renormalized masses m*p. The theoretical 
value for T\ is very close to the experimental result of 2.18 K, if the choice m*p ^ 1.6m 



is adopted. The results obtained for the connection between the predicted T\ and 
the effective quasiparticle mass are essentially independent of the value chosen for 
the mass parameter m*oii. We note further that a judiciously tailored temperature- 
dependent m*p(T) could be used to match the theoretical results for the specific heat 
to the experimental data in the normal phase. 



CONCLUSIONS 

In summary, we have performed a complete functional optimization of the den- 
sity matrix of liquid ^He, based on a trial form of generalized Jastrow type involv- 
ing only two-point descriptors of the dynamical correlations and the collective and 
single-particle thermal excitations. The Euler-Lagrange equations of this approxi- 
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mate description have been solved numerically in the interesting temperature range 
< T < 5.5 K. The ansatz for the density matrix admits a simultaneous treatment 
of the Bose-Einstein condensed phase and the normal phase. However, in the ab- 
sence of backflow effects, the predicted lambda transition temperature of Tx ~ 3.4 K 
exceeds the experimental value of 2.18 K by some 1.2 K. We argue that this short- 
coming of the theory stems from an inability of the chosen trial density matrix to 
accommodate important variations of effective masses - an inflexibility that results 
in high-lying collective excitations and quasihole and quasiparticle excitations that 
drive the theoretical transition temperature to unrealistically high values. This view 
is strengthened by the results obtained when the optimization is repeated in the nor- 
mal phase with an adjustable effective mass for quasiparticles. A renormalized mass 
m*p = 1.6m leads to a transition temperature nearly coincident with the experimen- 
tal value. 

To pursue a systematic microscopic inclusion of backflow and effective mass ef- 
fects, the present treatment should be generalized to allow for triplet pscudopotentials 
or three-point factors in the temperature-dependent wave function (2) and the inco- 
herence factor (3). These extensions will be required to achieve quantitative accuracy 
for the correlation energy and the dispersion relation of phonons and rotons. Existing 
techniques [44,45] may be appropriated to carry through the associated analysis and 
computations, which should include effects of elementary diagrams. In addition, a 
generalization of the ansatz (4) will be needed to take proper account of the quasipar- 
ticle and quasihole effective masses. The introduction of shadow coordinates [42,43] 
offers an alternate route to successful microscopic extension of the current realization 
of density-matrix theory. 
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In this appendix, we collect the necessary formulae for energies, occupation 
numbers, generalized structure functions, and related quantities. 

The ingredients £'c, ecoii, Ccc of the energy functional defined by Eqs. (5)-(8) 
have the explicit expressions 
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cc 



(k) l + Sif{k)+2B,,S^^\k) , (Al) 



where 



v*{r) = v{r) 



4m 



Au{r) 
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is the Feenberg effective potential [9,10], and 



ecoii(/i;) = le,{k)\ [l +47(/c)^(/c)l'^' + 1 \Sik)-^ 



(A2) 



ecc{k) = eo{k)\ 1 - f 1 - - Tcc{k)) kc(fc) - ^cc(fc) 



(A3) 



Tfie tfiermal occupation numbers associated witfi collective, quasihole, and quasipar- 
ticle excitations are given respectively by 



n{k) = ^ 



1/2 



(l + 4'y{k)S{k) 
n,,{k) = T,,{k)[l + Sif{k) 
n,{k) = V,,{k) \l + SfXk) + 2B,,S^''\k) 



(A4) 
(A5) 
(A6) 



where delivers F itself if F > and zero otherwise. The structure functions 
S{k) and Scc{k) and their various components, as well as the corresponding Fourier 
inverses g{r) = 1 + G{r) and Gcc{f), etc., are defined in Ref. [25]. They may be 
evaluated by means of the hypernetted-chain equations contained therein. 

The functions D{k) and Dcc{k) appearing in the Euler-Lagrange equations (15) 
and (16) are constructed as [28] 



dBr 



(A7) 



5rec(k) 



Iuj,,{k) + X] [l + SfXk) 



1 



+ -cc;,(fc)<(l + 5g)(/c) + 2S,, Sdd{k) + 2Sac{k) \Q{n,{k)) 



-eoik) 1 + Sacik) 
d 



5(°)(/c) + S, 



dBr 



Sdd{k) + 2Sdc{k) + SilHk)] }. (A8) 



The quantity Soik) = (Sddik), Sde{k), Seeik), Sddk), Sedk), Scdk)) is defined by 
the indicated sextuple of components. In Eqs. (A7) and (A8) these components form 
a scalar product with the six components T*^*) of the quantity T{k),i = 1,2, 3, 6. 
The component T^^^ is made up of two terms. 



T^^\k)=Ti^\k)+T^^\k) 



(1)/ 



(A9) 
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the first being the Fourier transform 



of 



VifiVi«)(ri2) + ree(ri2) 



+ ViiVed(ri2) (^4S2^iVde(ri2)ViiVed(ri2) 

+ 2Scc(l + iVdc(ri2)) [2ScciVdc(ri2)ViiVcd(ri2) 
+Vi(iVi2)(ri2) + r,e(ri2))" 
+ Vi [m\ri2) + rec(ri2)) [(l + iVdc(ri2)) Vi7Ved(ri2) 
+Vi(iVi?(ri2) + rec(ri2) 



and the second being given by 



T«(/c) = ^5ee6o(/c) 



2B,Ad{k) + SfXk)] \Xdc{k) + X,,{k) + Xf^k) + TcclA;) 



+ Sccrcc(/i;) tOcc{k) - i^c{k)Q{nc{k)) 
-]^B,Mk)[xfXk) + T,,{k)^ . 

For T(2) and T^^) we have 



1 - Xdd{k) - Xdeik) - 2B,,Xdc{k) 

- X,d{k) ^BecScd{k) + ^(^) Jxdc(A;) + Xec(/c) + X^^) (fc) + r 

- 'sfXk)Xdd{k) + (2Sec<5ed(A;) + SfXk))Xdc{k)\ 
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and 

(AM) 

The component T^^) again consists of two parts, 

T^''\k) = Ti''\k)+T^''\k), 
the first being the Fourier transform 



(A15) 
(A16) 



of 



X {25,,(l + iV,,(ri2)) ViiVed(ri2)+iVi2)(ri2)+r,,(ri2)} (A17) 
and the second having the expression 



SfXk) (l - Xdd{k) - Xde{k) - 2BccXdc{k)) 

- {2B,Ad{k) + SfXk)) {Xdcik) + X,,{k) + Xf^k) + T,,{k)) 
X {2B,,X,d{k) + X'g\k) + rec(A;)) 

2Sec^ed(/c) + ~SW{k)) (l - 2B,,Xdc{k) - Xil\k) - T,,{k) 



-2B,,~SfXk){x^d{k) + Xdc{k)) 



Xcd{k) + Xee{k) + Xjil\k) + Tec{k) 



+ 2B, 



2B,,S,d{k) + Sifik) 



Tcc{k){uJcc{k) - iOc{k)Q{nc{k)) 

-^eo(fc)(xW(fc)+ree(/c))' . 

The fifth and sixth components read 

T(^)(/c) = ^Seceo(A;)|xed(A;) 

X (l - 2B,,Xdc{k) - xfXk) - rcc(A;) 
-2B,,SfXk) (Xddik) + Xdcik)) 
- (sg\k)Xdd{k)+(2BccScd{k) + Sil\k))Xdc{k) 



(A18) 



X (2BccXcd{k) + Xg\k) + Tccik)) I , 



(A19) 
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X (l - 2B,,Xdc{k) - Xg\k) - r,e(A;)) 

-2B,J^^Xk) [Xdd{k) + Xac{k)) \ {2B,,X,d{k) + X^f{k) + r,,(A;)) 



1 

+ 2 



rcc(A;)(^Wcc(/f;) - a;c(/c)6(nc(/c))^ 
leo{k)[xi^\k) + TUk)y 



(A20) 



For the normal phase {Bcc = 0), these expressions simplify considerably, since many 
terms vanish [28]. 

The generalized structure functions S{k) and Scc{k) entering the Euler-Lagrange 
equations (15) and (16) may be explicated using the general prescription [28] 



+ Ti Eroo(*^')(i'io(*;') - i^cc(k') - a) 

k' I 

4.o(«=')(^<?('c')+r„(.'))1^, (A21) 



where a(3 = dd, de, ee, ec, dc, or cc (see Rcf. [25]). The dotted function Sa/3 may 
be analyzed and evaluated within the extended Fermi hypernetted-chain formalism 
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